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A new mechanism of tunnelling at macroscopic distances is proposed for a wave 
packet localized in one-dimensional disordered potential with mirror symmetry, 
V{—x) = V{x). Unlike quantum tunnelling through a regular potential barrier, 
which occurs only at the energies lower then the barrier height, the proposed mecha- 
nism of tunnelling exists even for weak white-noise- like scattering potentials. It also 
exists in classical circuits of resonant contours with random resonant frequencies. 
The latter property may be used as a new method of secure communication, which 
does not require coding and decoding of the transmitting signal. 

PACS numbers: 72.15.Rn, 72.10.Bg 72.80.Ng 

It is well-known that all quantum states in one- dimensional white-noise potential are 
strongly localized and quantum transport is limited by the distances not exceeding the 
localization length 1{E). At longer distances the destructive interference between direct and 
backscattered waves suppresses exponentially the amplitude of a wave packet. Statistical 
correlations in the disordered potential may change the interference pattern and may give 
rise to a discrete set [1] or to a continuum of delocalized states {2, 3] for short- or long-range 
correlation respectively. Correlations is a manifestation of the local properties of a random 
potential. The symmetry is a global property, therefore its effect on the transport may be 
even stronger. 

In this Letter we propose a symmetry-driven mechanism of tunnelling, which is specific 



for the random potentials only. Usually, the symmetry is considered to be an irrelevant prop- 
erty in disordered systems since the wave functions are localized at the (local) scales, which 
are much smaller than the (global) scales, where the symmetry is manifested. However, the 
symmetry of the potential, V{—x) = V{x), leads to definite parity of the wave functions. 
Either parity (even or odd) of an eigenfunction means that there are two equal peaks with 
half-width ~ 1{E) centered at the symmetric points. A symmetry- induced correlation be- 
tween these peaks gives rise to the mechanism of tunnelling of a wave packet (or excitation), 
independently how far apart the peaks are. Due to this mechanism a wave packet tunnels 
at macroscopic distances - a process which does not exist for the random potential without 
the symmetry. Natural disorder usually does not exhibit the mirror symmetry. Neverthe- 
less, the proposed mechanism of tunnelling is not of pure academic interest, since it may be 
observed also in a classical system - a random electrical circuit, where the symmetry can 
be easily introduced. In what follows we propose a new method of secure communications 
based on the symmetry-induced mechanism of tunnelling. The merit of this method is that 
it does not require a coding- decoding procedure. 

To demonstrate the main idea of the symmetry-induced tunnelling we consider the tight- 
binding Anderson model {J] . For one-dimensional lattice a stationary solution for the eigen- 
state with energy E is obtained from the equation 



where e„ is on-site energy. The energies E and e„ are measured in units of the hopping 
amplitude t, which in the case of diagonal disorder is independent on the site index n. 

Discrete Schrodinger equation ([T]) gives exact description of the electrical circuit of clas- 
sical impedances Z„ and z„ shown in Fig. [H Application of KirchhofT's Loop Rule to three 
successive unit cells of the circuit leads to the following linear relation between the currents 
circulating in the (n — l)-th, n-th and {n + l)-th cells 



If the vertical impedances are all the same, Zn = zq, this equation is reduced to the tight- 
binding model with diagonal disorder, Eq. ([T]), with e„ = 6n/zo and E = 2 + Zq/ zq. Here the 
random value of the impedance is split into its mean value Zq = and the fluctuating 
part 6n = Zn- Zq. 
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FIG. 1: Segment of infinite electric circuit of classical impedances. 



This exact correspondence allows testing of quantum effects of Anderson localization 
using classical electrical circuits with random elements. In fact, during the last decade 
chaotic resonant cavities have been successfully used for testing the predictions of quantum 
chaos j^. It is worth mentioning that electrical circuits of passive elements have been 
widely used for modelling different physical phenomena. The first application of the method 
of equivalent circuit probably goes back to Lord Kelvin who used a discrete RC chain to 
study a signal transmission through a transatlantic cable. Many bright examples of electrical 
circuits that model quantum mechanical behavior for simple but fundamental systems are 
g,ve„ the book by Pippard g Recently it was proposed that electro—tie waveguide 
can be used to model as exotic effect as Hawking black hole radiation [7i]. Some effects 
of correlated disorder have been studied in the experiments with microwave propagation 
;hrough disordered waveguides {s] and sub-terahertz response of superconducting multilayers 
91] . Experimental realization of a system with desirable correlations and observation of the 
localized and extended states are much easier in electromagnetic devices than in 

heterostructures with intentionally introduced disorder 10 1. 

If the potential in Eq. ([T]) is an even function, e„ = e_„, the eigenfunctions are either 
even or odd functions of n. If an eigenfunction is localized near a site no, the amplitude 
of this state at the origin is exponentially suppressed, '^n=o exp(— | uq \ /l{Eaj), provided 
I no 1^ l{Ea). However, due to definite parity of the wave function, another peak appears 
at the symmetric point n = —tiq. Strong localization of any excitation in random potential 
is a result of destructive interference between propagating and backscattered wave. The 
appearance of the symmetric peak can be explained as a result of constructive interference. 
It leads to exponential increase of the amplitude of the wave, i.e. to antilocalization 111 ]. 

In Fig. [2] we show two quasi-degenerate eigenstates calculated for the symmetric potential 
of 1000 sites (i.e. only 500 of these sites are random). The inverse localization length (the 




FIG. 2: (Color on line) Two eigenstates with different parity {^s is even and is odd) in random 
symmetric potential e_n = e„ with (e„) = and (e^) = eg = 0.1. These states belong to a doublet 
with energy splitting AE. Inserts: (a) Blow-up of the right peaks of the eigenfunctions showing 
that they possess different parity, (b) Numerical result for the localization length as compared to 
the energy- independent function 1{E) = 40. The compensation of the energy dependence in lo{E) 
is not of principal importance and is done only to simplify the discussion of the numerical results. 

Lyapunov exponent) can be estimated from the formula 31] 

oo 

l-\E) = l,\EMf^), ^(/x) = l + 2^eWcos(2/xfc). (3) 

A;=l 

Here Iq^{E) = eo/(8sin^/i) is the Thouless [3] result for the white noise disorder, 
the function v^(/i) accounts for the contribution of correlations with correlation function 
{en€n+k) = €-l^{k), and the dispersion relation is E = 2 cos /x. The results shown in Fig. [2] are 
obtained not for white-noise but for slightly correlated disorder with correlator ,^(1) = —1/2 
and C,{k > 1) = 0. These short-range correlations are introduced in order to compensate the 
smooth energy dependence of lo{E). It is easy to see that the contribution of the term with 
= 1 in Eq. ([3]) provides a flat dependence 1~^{E) = e^/A = const. Insert (b) in Fig. [2] 
shows the numerical values of 1{E), which fluctuate around 40 - the value obtained from Eq. 
([3]). In agrement with this estimate, the half- width of the peaks in Fig. [2] is approximately 
40 sites. 

The energy spectrum of Eq. ([T]) with symmetric random potential is similar to the 
spectrum of a double-well potential. It consists of discrete levels, most of them lying within 
the interval —2 < E < 2. The energy levels are arranged in doublets of states with different 
parity. The energy 6{E) between the centers of the doublets scales with the length of 



the system as The energy sphtting /S.E in the doublet is exponentially small, 

A(£^) oc exp[— 4 | uq \ //(-E)], i.e. the states are quasi-degenerate. Both, 5{E) and A(i?) 
fluctuate with energy because of statistical fluctuations of the density of states, uq^E), and 
1{E). 

The symmetry-induced tunnelling can be observed in the dynamics of an excitation. Let 
a perturbation is applied at one of the sites of the symmetric random sequence. In the 
simplest case the perturbation is a 5-excitation at the site uq, ipnit = 0) = 5„„p. Since this 
excitation is not an eigenfunction of the system, its temporal evolution is represented as a 
superposition, 

Mt) = J2c:o'^n^M-^E^t)- (4) 

a 

The sum in Eq. runs over the eigenstates, which are all localized. The eigen- 

states centered closer to the initial excitation contribute more because the coefficient 
^no ~ {'^nl'^nit = 0)) = "^^^ is the Overlapping integral between the initial excitation 
and the eigenstate Let the eigenstates with maximum overlapping be and "^s- They 
form a doublet with the central energy E = {Ea + Es)/2 and splitting AE = Ea — Es- 
Taking into account only these two terms in Eq. (jlj), the following approximate result for 
the evolution of the initial excitation can be easily obtained 

-iEt r //^^^ X ... /AE 



COS I -^t I + i sin [^-^t I ^-(^1 



(5) 



Here \l/±(ri) = ("^ a) / V^- Each of these linear combinations is a single-peak function. 
The peak of \E'4. is always close to the point of initial excitation. For the eigenf unctions 
shown in Fig. [2] the peak is localized in the region of negative n. 

At the early stage of evolution the initial 5-peak at uq quickly spreads over the region of 
width 21{E). Further spreading is suppressed by Anderson localization. Eq. is not valid 
at this transient stage. This equation describes steady and "slow" harmonic oscillations of 
the initial excitation between the two symmetrical points. The period of oscillations of the 
density | ipnit) P is T = 27t/AE. If the distance 2 | no | between the peaks exceeds the 
localization length, the amplitude of the wave function at the origin is exponentially small 
~ exp(— I uq I /1{E)). It, however, grows exponentially towards the symmetrical point —uq. 
This increase is a manifestation of the tunnelling induced by the symmetry. The dynamics of 
penetration of the initial excitation to the symmetrical point is very similar to the tunnelling 
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FIG. 3: (Color on line) Temporal evolution of the probability density | tp{t) ^ calculated at the 
point of excitation, no = —153 (grey line) and at the symmetric point, —uq = 153 (black line). 
Dashed lines are the window-average values of the densities at the symmetric points and are drawn 
for guide-eye only. The insert shows similar temporal evolution but for much stronger disorder, 
eg = 0.55. In this case the localization length is shorter 1{E) = 8 and much less doublets contribute 
to the evolution of the the initial excitation. 

through a potential barrier, although there is no real barrier. Exponential decrease (increase) 
of the wave functions is due to multiple scattering events with predominant destructive 
(constructive) interference. One can speak about an effective double-well potential which 
produces the same discrete energy spectrum. Calculation of the parameters of this effective 
potential is a challenging inverse-scattering problem. Tunnelling processes without a real 
barrier are known in dynamical systems, where quantum transitions occur either between 



strongly localized states [13|| or between classically separated regions in phase space [M]- It is 
worth mentioning that regular Bloch-like oscillations in a potential with correlated disorder 
may occur also due to the presence of two mobility edges in the energy spectrum [l^. 

Eq. takes into account interaction of the initial excitation with the nearest doublet. 
If there are more eigenstates in Eq. (jl]), whose wave functions extend to the point no, they 
also contribute to the evolution of the initial excitation. In this case the oscillations between 
the peaks at uq and —uq are not harmonic any more but a superposition of harmonics with 
different periods. In the numerical study of evolution of the excitation we take into account 
its interaction with the eigenf unctions which have amplitude > 10^'^ at the site uq = —153. 
There are 540 such eigenfunctions out of total 1000. These eigenf unctions produce the 
oscillatory pattern in Fig. El Although the site Uq = —153 is the position of the maximum 
for the eigenstates and in Fig- El other states give a noticeable contribution. Since 
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FIG. 4: (Color on line) Spatial distribution of the probability | Vn(i) P at three instants: (a) 
Spreading of the initial peak at the transient stage, t = 9.6- lO^^T. The initial peak of amplitude 1 
at no = —153 is shown in grey. The secondary peak is not visible at this stage, (b) The secondary 
peak at —no = 153 is well developed at t = 0.241T; (c) The two peaks become almost equal at 
t = 0.607r. 

the levels splittings in different doublets are random and incommensurate, the evolution of 
the wave packet is not periodic but it keeps the main features predicted by Eq. In 
the case of stronger localization of the eigenstates the dependence | ^Jnoit) p approaches the 
harmonic dependence Eq.Q as it is seen in the insert in Fig. El 

Spreading and tunnelling of the initial excitation is shown in Fig. |H At the transient 
stage the excitation broadens up to the size of 21{E) ^ 80, Fig. Hti. The initial stage is 
followed by the long-lasting stage of tunnelling at the macroscopic distance 2 | n-o |. The 
tunnelling gives rise to the secondary peak at —no, which "slowly" grows and reaches its 
maximum at t ^ Fig. Hl^. The amplitude of the secondary peak in Figs. [S] and 
does not exceed 5% of the initial peak, but at t ^ T/2 the both peaks have approximately 
the same amplitude. The amplitude of the peaks can be obtained from the normalization 
condition and it is determined by 1{E) as it follows from Eq. ([5]). This amplitude is much 




larger then the amphtude at the origin n = 0, as it can be clearly seen from Fig. |H The 
amplitude of the peaks at the points ±nQ increases with eo- Simultaneously the wave function 
at the origin decreases exponentially and can be easily controlled by the disorder. 

Application of the proposed ideas to random electrical circuits is straightforward. In 
what follows we demonstrate the evolution of the signal in a circuit shown in Fig. [1] with 
vertical impedances being equal solenoids with and horizontal impedances 

being capacitors with Z„ = ^ (^1 — Here 5C„ is the fluctuating part of the 

capacitance, which is an even function of n, 6Cn = SC^n- Propagation of an excitation in 
this circuit follows the wave equation 

Cn+lVn+l + Cn-lVn-l — 2CnVn = Vnj -^0; (6) 

where Vn is voltage drop at the n-th capacitor. This equation requires two initial conditions. 
Let the voltage drop Vq is applied at t = to the capacitor C„o, inducing the initial current 
-^0 = CnoK)- Stationary solutions (V^ oc exp(— icut)) are either even or odd functions and the 
spectrum of eigenfrequencies consists of a set of doublets. For an infinite chain the majority 
of the eigenfrequencies occupy an interval [cuq/S, oo], where = [Li^C^ ' . Assuming that 
the initial perturbation excites only the closest to the site pair of eigenstates {Va and 
Vs)^ the solution of Eq. (I6l) can be written in the form similar to Eq. ([5]), 
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C„o Vo sin(co't) — — cos(co't) 



C„,)Vo cos (cut) + — sin(ct;)f:) 



sin (^t)^-H 



cos(^— (7) 



Here uj is the center of the doublet, and Acj is the frequency splitting. Single-peak functions 
^±(^) = (^5 ± Va)/\/^ play the same role as do in Eq. ([5]). Eq. ([7]) shows that 
the evolution of the initial signal in a random (symmetric) electrical circuit is similar to 
the wave packet evolution obtained from the tight-binding model. There is an obvious 
symmetry- induced tunnelling of the initial signal at macroscopic distances. 

The effect of tunnelling can be used for secure communications. Instead of coding and 
decoding a signal, we propose to suppress the transmitted signal by a circuit with random 
elements and then to restore it, using a symmetric counterpart of the random circuit. The 
signal can be suppressed to the noise level and safely transmitted to the receiver over a 
transmitting line. The symmetric counterpart of the random circuit restores only the signal, 



(not noise) since constructive interference occurs only for the coherent part, which has 
passed trough the suppressing circuit of the emitter. The non-coherent part (noise or any 
irrelevant signal) will be exponentially suppressed by the receiving random circuit. The two 
identical random circuits may be fabricated as microchips, which are installed (or replaced) 
simultaneously at the emitter and receiver. In the absence of dissipation and asymmetry 
between the two random elements, the proposed method guarantees a robust restoration of 
the signal. Inevitable Joule losses should be compensated by an amplifier, which does not 
destroy the coherency of the signal. 

In conclusion, we demonstrate that in a symmetric random potential the localized quan- 
tum states have two peaks as it is required by the parity. Fast spreading of the initial 
(5-excitation within localization length is followed by slow growth due to tunnelling at the 
symmetrical point. This effect opens a new possibility for secure communications that does 
not require coding and decoding of the transmitting signal. The random circuits may oper- 
ate in a wide range of radio-frequencies, using commercial capacitors and inductors. They 
also can be fabricated and operated in the infrared and optical region using the concept of 
plasmonic nanoelements proposed in Ref. 16 1. 
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